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1.  INTRODUCTION 


Almost  all  flows  of  interest  to  ballisticians  involve  the  presence  of  a  second  phase  and 
turbulence.  In  the  absence  of  theoretical  or  experimental  insight  into  their  dynamics,  the 
adopted  approach  has  been  either  to  ignore  the  turbulence  or  use  ad*hoc  single-phase 
descriptions  with  a  number  of  "adjustable”  parameters.  Many  of  these  do  not  have  a  sound 
theoretical  underpinning  or  involve  extrapolations  and  are  more  In  the  nature  of  post-facto 
instead  of  predictive. 

To  establish  better  modeling  and  computational  capability  in  this  field,  a  research  effort 
was  undertaken  to  shed  light  on  the  interaction  of  par ':culate  presence  on  turbulence  and  vice 
versa.  The  statistical  theory  of  turbulence  gives  a  reasonably  accurate  descriptiori  of  the 
measured  average  flows  and  correlation  functions  and  can  be  extended  to  nonreactive  two- 
phase  flows.  However,  if  the  two  phases  consist  of  a  gas  and  of  a  dispersed  solid  phase,  the 
extension  is  not  a  straightfonvard  one.  In  fact,  it  requires  a  certain  amount  of  "coarse- 
graining."  (i.e.,  spatial  averaging);  this  question  was  discussed  rather  carefully  by  Besnard  and 
Harlow  (1988).  Alternatively,  one  can  describe  tfie  dispersed  particulate  phase  by  means  of  a 
Boltzmann  equation  (Domokos,  Kovesi-Domokos,  and  Zoltani  1988a,  1988b).  This  has  the 
advantage  that,  in  terms  of  a  Chapman-Enskog  expansion,  one  can  generate  not  only  the 
Eulerian  equations,  (the  result  of  coarse-graining)  but,  in  principle,  corrections  of  arbitrary 
order  to  it.  The  difficulty  is  that  very  soon  one  runs  into  substantial  computational  difficulties. 
Even  with  the  help  of  modern  computers,  a  straightfonA/ard  approach  to  solving  the  resulting 
coupled  set  of  equations  can  be  very  time  consuming. 

In  a  previous  work  (Domokos,  Kovesi-Domokos.  and  Zoltani  1991)  we  proposed  a  series 
of  variational  principles  in  the  framework  of  the  statistical  theory  of  single-phase  turbulent 
flows,  based  upon  the  work  of  Martin.  Siggia,  and  Rose  (1973)  and  De  Dominicis  and  Peliti 
(1978)  (see  also  Domokos,  Kovesi-Domokos,  and  Zoltani  1988a.  1988b).  A  variational 
principle  has  the  advantage  that  there  exist  methods  to  extremize  the  functional  in  question 
which  are  very  economical  from  the  calculational  point  of  view— the  Rayleigh-Ritz  method 
being  the  most  notable  one.  (The  disadvantage  is.  of  course,  that  some  insight  is  necessary 
in  order  to  guess  good  trial  functions;  most  of  the  time  one  proceeds  by  trial  and  error  in 
several  steps.)  In  Domokos,  Kovesi-Domokos,  and  Zoltani  (1991)  we  performed  a  sample 
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alculation  in  order  to  describe  a  singie-phase,  cyiindrically  symmetric  jet.  Despite  the 
implicity  of  the  triai  functions  and  the  smali  number  of  parameters,  quite  a  reasonable 
igreement  with  the  data  was  achieved. 

The  purpose  of  the  present  work  is  to  extend  the  calculation  to  nonreactive,  two-phase 
lows.  The  plan  of  this  report  is  the  following:  In  the  next  section  we  briefly  review  the 
ormalism  described  in  Domokos,  Kovesi-Domokos,  and  Zoltani  (1988a,  1988b,  1991).  For 
ho  sake  of  brevity,  we  use  a  somewhat  abstract  notation — ^this  makes  the  concepts  and  the 
structure  of  the  formalism  more  transparent.  In  Section  3  we  state  the  framework  for  the 
somputation  of  two-phase  flows.  Section  4  contains  a  sample  of  calculations  of  some  of  the 
jroperties  of  a  cyiindrically  symmetric  two-phase  jet  and  a  comparison  of  the  results  with 
3xperimental  data.  Finally,  Section  5  contains  the  conclusions. 

2.  THE  FORMALISM 

We  consider  a  vector  space  of  dynamical  variables;  X  denotes  a  generic  element  of  the 
t/ector  space.  In  general,  X  is  a  function  of  space  and  time.  For  instance,  X  may  stand  for 
the  six  components  of  the  velocity  field  of  a  two-phase  flow  at  a  given  space-time  point  (i.e., 
three  components  of  the  velocity  of  the  carrier  fluid  and  three  components  of  the  velocity  of 
the  particulate  component).  We  assume  that  X  obeys  an  autonomous  equation  of  motion  of 
the  form 


d,X^F[X]  =  f,  (1) 

where  the  functional  F[X]  may  contain  spatial  derivatives,  integrals  over  spatial  coordinates, 
Btc.,  but  no  integral  over  time  and  no  time  dependence  either.  Further,  f  is  a  Gaussian 
random  force,  with  correlation  operator  K. 

The  generating  functional  of  the  correlation  functions  is  given  in  terms  of  an  arbitrary 
source,  j ; 


m-  fDxexp[-<d,X  -  F\K{d,X  -  F)>  *  <j\X>]  .  (2) 
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Here  <  .1  .>  stands  for  a  suitable  scalar  ptfoduct  over  the  vector  space,  including  integration 

I  '  .  . 

over  space-time  variables  and  summation  over  discrete  components.  The  functional  measure, 
jOx.  contains  an  infinite  determinant,  viz.  det  ^d,  -  5F/SX;,  as  expistined  in  Domokos,  Kovesi- 
Domokos,  and  Zoltani  (1988a).  The  cumulants  are  generated  by  IV  =  -In  Z. 

The  various  averages  are  obtained  by  taking  functional  derivatives  of  W  with  respect  to  j. 
We  use  the  notation 


G(1)h<X(x,)>  = 


&W 

Sy(x,)  ’ 


(3) 


G(1.2)  =  <X{x,)X{x^)> 


5‘W 

^j{x,)Sj{x^)  ’ 


(4) 


and  so  on.  Here  Xj  stands  for  a  space-time  point. 

in  order  to  get  a  suitable  variational  principle,  we  also  add  a  bilinear  source  in  the 
exponential  of  Equation  2,  of  the  form,  <  q  (1,2)  I  G(2, 1)  >  (Domokos,  Kovesi-Domokos,  and 
Zoltani  1991).  Next,  we  perform  a  double  Legendre  transformation  in  the  variables  y  and  t|,  so 
that  the  resulting  functional  has  G(1)  and  G(2)  as  its  functional  arguments.  We  denote  this 
functional  by  S. 

One  has  the  relations 


5S 

8G(1) 


-y(i) 


5S 

5G(1,2) 


-ti(1.2)  . 


(5) 


The  functional  given  by  Equation  5  is  stationary  if  the  arbitrary  sources  are  put  equal  to  zero. 
The  reader  will  readily  recognize  that  functionals  of  the  type  in  Equation  5  play  a  role 
analogous  to  the  entropy  in  statistical  mechanics.  For  this  reason,  relations  of  the  type  5  with 
vanishing  external  sources  were  called  the  principle  of  stationary  entropy  in  an  analogous 
context  by  Oe  Dominicis  and  Martin  (1964). 
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3.  EQUATIONS  FOR  A  TWO-PHASE  TURBULENT  FLOW 


In  what  follows,  variables  characterizing  the  carrier  fluid  w\\\  be  given  a  subscript  f,  those 
characterizing  the  dispersed  particulate  phase  a  subscript  p.  We  work  in  terms  of 
dimensionless  variables  by  dividing  velocities  with  some  characteristic  speed,  coordinates  by  a 
characteristic  size  of  the  system  under  consideration,  etc.  In  this  report  we  write  down  the 
equations  of  motion  in  the  leading  approximation  of  the  Chapman-Enskog  expansion 
(Domokos,  Kovesi-Domokos,  and  Zoitani  1988a}  so  that  both  phases  obey  the  equations  of 
hydrodynamics.  The  volume  fractions  of  the  fluid  and  particulate  phase  are  denoted  by  e^and 
Ep.  respectively.  Of  course. 


Zf  +  tp=1  . 


(6) 


We  now  have  the  equations  of  continuity. 


a,£, .  ^U,ut]  -  0 ,  (7) 

oX 

S.e.  *  -  0  •  (8) 


The  equations  of  motion  read 


Here  /}and  fp  stand  for  the  perturbing  Gaussian  random  forces  acting  on  the  fluid  and 
particulate  phases,  respectively.  In  Equations  9  and  10,  p  stands  for  the  external  pressure. 
The  quantity  C,  in  general,  is  a  function  of  I  Uf  Up].  However,  with  the  exception  of  some 


unusual  cases  such  as  highly  viscous  caniers,  very  heavy  loading,  etc.,  the  velodty  difference 
between  both  phases  is  not  too  big.  In  that  case,  C  can  be,  to  a  good  approximation, 
replaced  by  its  value  given  by  Stokes’  taw.  This  results  in  a  considerable  simplification  of  the 
computations. 

The  reader  will  notice  that  a  viscous  term  has  been  omitted  from  Equation  9.  Despite  the 
fact  C  is  proportional  to  the  viscosity  of  the  carrier  fluid,  the  approximation  is  a  permissible 
one  unless  one  is  interested  in  very  small  scales  (large  wave  numbers).  The  coupling  temi  is 
proportional  to  a  velocity  difference  (in  the  Stokesian  approximation),  whereas  the  viscous 
term  is  proportional  to  the  scalar  curvature  of  a  velocity  field  (on  large  and  moderate  length 
scales,  the  latter  is  less  important  than  the  former). 

Let  us  now  take  a  look  at  Equations  6.  7.  and  8.  One  immediately  realizes  that  if  the  flow 
of  one  of  the  phases  is  approximated  by  an  incompressible  one,  the  flow  of  the  other  phase 
becomes  incompressible  too.  In  a  large  number  of  practically  important  situations,  the 
approximation  of  an  incompressible  flow  is  a  rather  good  one.  In  what  follows,  we  are  going 
to  make  the  approximation.  Approximating  a  flow  by  an  incompressible  one  has  two 
immediate  consequences.  First,  the  pressure  is  no  longer  an  independent  dynamical  variable; 
it  can  be  integrated  out  explicitly  from  the  generating  functional  of  the  correlation  functions.* 
Second,  the  equations  of  continuity  now  tell  us  =  V-u^  =  0;  hence,  the  velocity  fields  can 
be  obtained  as  curls  of  vector  potentials  with  an  ensuing  gauge  of  freedom  (Domokos,  Kovesi- 
Domokos,  and  Zoltani  1991). 

Finally,  the  expression  of  the  stirring  force  (or,  more  precisely,  its  correlation  operator)  has 
to  be  discussed.  In  Section  4,  and  in  many  other  practically  important  cases,  we  are 
concerned  with  problems  where  the  mean  flow  is  cylindrically  symmetrical,  with  the  mean 
velocity  having  a  large  component  along  the  axis  of  symmetry  and  a  rather  small  radial 
component.  In  such  cases  one  gets  satisfactory  results  by  taking  a  stirring  force  of  the  same 
symmetry  and,  in  fact,  neglecting  the  radial  component  of  the  stirring  force  altogether.  We 
take  for  both  phases  an  identical  form  of  the  matrix  elements  of  the  correlation  operator,  viz., 


*  Ttad'rtbnally,  this  is  formulated  as  using  the  equations  of  motion  to  eliminate  the  pressure;  however,  within 
the  present  context,  performing  a  Gaussian  functional  integral  over  the  pressure  leads  to  the  desired  result  in  a 
more  transparent  manner  (Domokos,  Kovesi-Domokos,  and  Zoltani  1991). 
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K,,,(x,)f')“exp(-alx»  -  x'»  |)5»(x,  -  x,'  )5(t  -  t')S,,,S,,,. 


(11) 


where  the  axis  of  symmetry  has  been  chosen  as  the  3^  axis  and  x^- denotes  a  vector  lying  in  a 
plane  perpendicular  to  it.  The  stream-wise  correlation  of  the  stirring  force  is  governed  by  the 
parameter  a.  As  it  turns  out,  Its  magnitude  is  not  a  critical  one  (Burgett  1989),  but  it  seems 
that  letting  a  3  is  not  a  very  good  approximation. 

Now  we  have  assembled  the  elements  of  computing  the  generalized  entropy.  Equation  5. 
We  have  done  it  to  two-loop  accuracy  (Domokos,  Kovesi-Domokos,  and  Zoltani  1991)  (i.e.,  by 
computing  the  first  approximation  to  the  solution  of  the  Cornwall- Jackiw-Tomboulis  functional 
differential  equation). 

Unfortunately,  the  result  of  the  calculation  is  neither  transparent  nor  revealing.  (In  fact, 
most  of  the  calculation  has  to  be  performed  by  means  of  symbolic  manipulation  programs.) 
The  reader  will  be  spared  the  sight  of  the  result.  Instead,  in  the  next  section,  we  present  the 
results  of  a  calculation  of  the  correlation  functions  of  a  steady  two-phase  jet  sufficiently  far 
from  the  plane  of  injection.  (In  this  way,  one  can  assume  that  the  turbulence  is  fully 
developed:  transients  died  away  and  the  correlation  functions  are  stationary.) 

4.  CORRELATION  FUNCTIONS  OF  A  CYLINDRICALLY  SYMMETRIC  TWO-PHASE  JET 

The  calculation  of  a  cylindrically  symmetric  jet  is  of  considerable  practical  importance — 
many  jets  possess,  to  a  very  good  approximation,  axial  symmetry.  In  addition,  it  is  a  relatively 
simple  configuration;  it  is  eminently  suitable  for  testing  a  method  of  calculation.  We 
proceeded  in  a  way  which  proved  to  be  successful  in  our  previous  work.  The  procedure 
follows: 

•  The  velocity  fields  are  obtained  as  the  curl  of  vector  potentials:  u,p  =  V  x  A^ 

We  work  in  an  axial  gauge:  the  component  of  the  vector  potential  along  the  axis  of 
symmetry  (chosen  to  be  the  3^  axis)  vanishes. 

•  A  Reynolds  decomposition  is  used  for  the  vector  potentials:  A  =  (A)  +  A' 
correspondingly,  tiie  velocity  fields  are  decomposed  as  u  -  U+tf.  (For  the  sake  of 
simplicity,  we  omitted  the  indices  f.p,  the  decomposition  is  used  for  both  phases.) 
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We  determine  the  form  of  the  correction  functions  of  the  vector  potentials: 


u 


(12) 


in  this  equation,  5  and  e  stand  for  the  Kronecker  and  Levi-Civita  tensors, 
respectiveiy.  Lower  case  Greek  indices  refer  to  vector  components  in  the  piane 
transverse  to  the  symmetry  axis,  the  superscripts  /,  j  assume  the  values  f  and  p. 

•  One  assumes  a  trial  form  of  the  vector  potentials  and  of  the  functions 
containing  a  few  unknown  parameters. 

•  These  expressions  are  substituted  into  the  expression  of  the  effective  Equation  5. 

•  The  entropy  is  then  extremized  witti  respect  to  the  parameters.  This  determines 
their  optimal  values  given  the  functional  form  of  the  vector  potentials  and  correlation 
functions. 

•  Finally,  the  mean  velocities  and  correlation  functions  are  computed  by  taking  the 
appropriate  curls. 


(Extremizing  the  entropy  with  respect  to  parameters  in  a  given  functional  form  is  the  Rayleigh- 
Ritz  method  of  solving  a  variational  problem.) 

All  velocities  are  measured  in  units  of  the  fluid  velocity,  Ug,  on  the  axis  at  the  plane  of 
injection.  The  unit  of  length  is  the  diameter  of  the  injection  pipe.  All  vector  potentials  and 
correlation  functions  are  assumed  to  be  time  independent.  We  chose  essentially  the  same 
trial  functions  for  both  phases  as  in  the  case  of  a  single-phase  flow.  (Of  course,  the  optimal 
values  of  the  parameters  are  different.)  We  used  the  3^  component  of  the  mean  velocity  on 
the  axis  as  an  input.  The  data  were  taken  from  Zoltani  and  Bicen  (1990).  The  following  forms 
give  a  good  fit  to  the  data. 


uUz.r^O)  » 
Ui{z,r«0)  -  0.78 


1.35 

1  +  0.035 z®'*  ’ 

(  1.13 

J  +  0.007 z®'* 


0.35  exp( -0.1  z®)  , 

-  0.13exp(-0.03z®)' 

> 


(13) 

(14) 
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The  functional  form  of  the  average  of  the  vector  potential  was  assumed  to  have  a  Gaussian 
radial  dependence  for  both  phases.  In  cylindrical  coordinates  we  have 


/ 

<A^)  -  l«(2)®exp 

V 


(15) 


with  R(z)  ^a  +  bz.  The  optimized  values  of  the  parameters  a  and  b  are 

a,  »  0.25,  b,  »  0.076  ;ap  =  1.38,  b^  =  0.013  . 

Likewise,  we  used  the  same  functional  form  for  tiie  quantity  Z,  as  in  Domokos,  Kovesi- 
Domokos,  and  Zoltani  (1991)  for  both  phases 


Z,  =  >texp[-/(z,  -  Zj)*  -  g{Xr,  - 

exp(-5M)(1  *  BM)[U^(z„0)U^{z^.0)y'’‘  .  (16) 


Here, 


2  2 
*  yi 

n(z.)“ 


+ 


Xj*  + 

R{z,y 


(Just  as  in  Domokos.  Kovesi-Domokos,  and  Zoltani  [1991],  we  set  Z^  =  0.)  The  optimal 
values  of  the  parameters  are 

A''  0.00^9,AP•>•  «  0.0013,8''=  2.3,8  "•»’  »  3.5,5''  =  1.79,5'''*  =  2.63  . 

The  values  of  the  parameters  f,  g  are  practically  identical  in  both  phases,  f'^-^  =  f''*~  5.6  and 
gp-Pm  g'-fm  3,7.  These  values  have  been  computed  for  the  particle  loading  extracted  from 
Zoltani  and  Bicen  (1990).  The  cross  correlation  functions.  Zf'^  have  basically  the  same 
shape  (the  parameters  8,5.f,g  are  practically  identical).  However,  the  overall  scale  is  much 
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smaller  «  A*’).  This  is  intuitively  envious— with  the  loading  fraction  used  in  Zoltani  and 
Bicen  (1990),  the  coupiing  between  both  phases  is  a  rather  weak  one.  As  a  consequence,  it 
is  rather  hard  to  obtain  a  reliable  estimate  of  AP‘*.  Several  iterations  fluctuate  around  a  value 
of  the  cross  correlation  scale  at  least  an  order  of  magnitude  smaller  than  (the  comparable) 
fluid  and  particle  correlation  scales;  however,  a  stable  maximum  of  the  generalized  entropy, 
Equation  5  is  hard  to  achieve. 

Once  we  have  these  parameters  values,  one  can  take  the  curl  of  the  mean  vector 
potential  and  the  double  curl  (with  respect  to  both  arguments)  of  ^  in  order  to  obtain  the 
mean  velocity  and  the  correlation  functions  of  the  fluctuations,  respectively.  Experimental  data 
on  the  fluctuation  correlation  are  often  taken  at  coincident  arguments;  this  was  the  case  in 
Zoltani  and  Bicen  (1990)  too.  We  display  the  results  for  those  correlation  functions  where 
data  were  available  in  those  articles.  The  other  correlation  functions  can  be  easily  reproduced 
from  the  preceding  formulae. 

In  the  figures,  we  return  to  a  conventional  notation  which  is  acceptable  for  coincident 
arguments.  The  correspondence  between  the  present  notation  (more  suitable  for  theoretical 
calculations)  and  the  conventional  one  is  the  following.  Define 


G...[^xZxVl..,  (17) 

(i.e.,  the  double  curl  of  the  fluctuation  correlation  function  of  the  vector  potential  with  respect 
to  both  arguments).  (In  the  case  of  coincident  arguments,  the  curls  are  taken  before  the 
a^uments  are  let  to  coincide.)  Note  that  defined  by  Equation  17  is  just  the  correlation 
function  of  velocity  fluctuations.  For  the  sake  of  simplicity,  we  omitted  the  indices  refer  to  the 
phase  {f,p)  in  question. 

We  now  have  the  correspondence  between  the  conventional  notation  and  ttie  one  used 
here  in  Domokos.  Kovesi-Domokos,  and  Zoltani  (1991) 
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(18) 


ui 


and  a  similar  relation  for  the  particle-particle  correlation  function.  In  an  arbitrary  Cartesian 
coordinate  system  (the  3^  axis  coinciding  with  the  axis  of  symmetry  and  the  orientation  of  the 

and  2^^  axes  being  determined  by  the  measuring  apparatus)  the  components  in  the  plane 
perpendicular  to  tfie  axis  of  symmetry  are  denoted  by  v  and  w,  respectively.  There  is  a 
relation  between  <v^.  <w^,  etc.,  and  our  notation  similar  to  the  one  exhibited  in  Equation  18. 

It  is  an  obvious  one  and  we  do  not  exhibit  it  here.  Note  that  cylindrical  symmetry  entails 

=  G2:4.Gi':4  =  o  .  (i9) 

for  any  combination  of  the  phases  p.f.  (There  is  no  good  notation  available  in  the 
conventional  system  of  notations  for  correlation  functions  such  as  Gf .  However,  among  the 
figures  presented  in  the  present  work,  no  such  correlation  function  has  been  plotted.) 

5.  CONCLUSIONS 

The  extension  of  the  calculations  described  in  Domokos,  Kovesi-Domokos,  and  Zoltani 
(1991)  to  a  two-phase  flow  have  basically  the  same  merits  as  in  the  case  of  a  single-phase 
flow.  One  notes  again  that  the  variational  method  is  a  very  economical  one  from  the 
computational  point  of  view.  With  rather  simple  functional  forms  of  the  trial  functions  and  in 
terms  of  a  few  parameters,  a  reasonable  agreement  with  the  experimental  data  can  be 
obtained  at  the  cost  of  a  relatively  small  computational  effort. 

One  can  conceive  a  number  of  ways  in  which  the  agreement  with  the  experimental  data 
could  be  improved. 

•  One  can  contemplate  calculating  higher  order  approximations  to  the  generalized 
entropy  (De  Dominicis  and  Peliti  1978). 

•  One  can  invent  (with  some  physical  insight)  more  sophisticated  trial  functions  for  the 
description  of  tfie  mean  flows  and  the  fluctuation  correlation  functions. 
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•  Clearfy,  more  experimental  data  are  necessary  In  order  to  determine  the 
dependence  of  the  various  correlation  hinctions  on  such  physical  quantities  as  the 
loading  fraction,  the  viscosity  of  the  carrier  fluid,  etc.  In  turn,  this  will  enhance  one's 
physical  insight  in  "guessing"  more  sophisticated  trial  functions  in  the  Rayieigh-Ritz 
method. 


At  this  point  it  is  hard  to  determine  the  increase  in  computational  cost  once  one  decides  to  go 
beyond  the  present,  simplest  approach.  Nevertheless,  experience  with  variational  methods  in 
various  branches  of  physics  suggests  that  many  methods  are  likely  to  be  rather  economical  in 
the  statistical  theory  of  turbulence  too. 
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Everett.  MA  02149-5936 

1  General  Electric  Company 

Tactical  System  Department 
ATTN:  J.  Mandzy 
100  Plastics  Ave. 

Pittsfield.  MA  01201-3698 

1  IITRI 

ATTN:  MJ.  Klein 
10  W.  35th  Street 
Chicago.  IL  60616-3799 

1  Hercules.  Inc. 

Allegheny  Ballistics  Laboratory 
ATTN:  William  B.  Walkup 
P.O.  Box  210 

Rocket  Center.  WV  26726 

1  Hercules.  Inc. 

Radford  Army  Ammunition  Plant 
ATTN:  E.  Hibshman 
Radford.  VA  24141-0299 

1  Hercules.  Inc. 

Hercules  Plaze 
ATTN:  B.M.  Riggleman 
Wilmington,  DE  19894 

3  Lawrence  Livermore  National 

Laboratory 
ATTN:  L-355, 

A.  Buckingham 
M.  Finger 

L-324,  M.  Constantino 
P.O.  Box  808 

Livermore.  CA  94550-0622 

1  Oiin  Corporation 

Badger  Army  Ammunition  Plant 
ATTN:  F.E.Wolf 
Baraboo.WI  53913 


No.  of 
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3  Oiin  Corporation 

ATTN:  E.J.  Kirschke 
A.F.  Gonzalez 
D.W.  Worthington 
P.O.  Box  222 

St.  Marks.  FL  32355-0222 

1  Oiin  Ordinance 

ATTN:  H.A.  McElroy 
10101  9th  Street,  North 
St.  Petersburg,  FL  33716 

1  Paul  Gough  Associates,  Inc. 

ATTN:  P.S.  Gough 
1048  South  St. 

Portsmouth,  NH  03801-5423 

1  Physics  International  Company 

ATTN:  LIbrary/H.  Wayne  Wampler 
2700  Merced  Street 
San  Leandro,  CA  984577-5602 

1  Princeton  Combustion  Research 

Laboratory,  Inc. 

ATTN:  M.  Summerfield 

475  US  Highway  One 

Monmouth  Junction,  NJ  08852-9650 

2  Rockwell  International 
Rocketdyne  Division 
ATTN:  BA08, 

J.  Flanagan 
J.  Gray 

6633  Canoga  Avenue 
Canoga  Park,  CA  91303-2703 

1  Thiokol  Corporation 

Huntsville  Division 
ATTN:  Tech  Library 
Huntsville,  AL  35807 

1  Sverdrup  Technology,  Inc. 

ATTN:  Dr.  John  Deur 
2001  Aerospace  Parkway 
Brook  Park.  OH  44142 

2  Thiokol  Corporation 
EHdon  Division 
ATTN:  R.  Biddle 

Tech  Library 
P.O.  Box  241 
Elkton,  MD  21921-0241 
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1  Veritay  Technology,  Inc. 

ATTN:  E.  Fisher 
4845  Mitlersport  Hwy. 

East  Amherst,  NY  14501-0305 

1  Universal  Propulsion  Company 

ATTN:  H.J.  McSpadden 
25401  North  Central  Ave. 

Phoenix,  AZ  85027-7837 

1  Battelle 

ATTN:  TACTEC  Ubrary,  J.N.  Huggins 
505  King  Avenue 
Columbus,  OH  43201-2693 

1  Brigham  Young  University 

Department  of  Chemical  Engineering 
ATTN:  M.  Beckstead 
Provo,  UT  84601 

1  California  Institute  of  Technology 

204  Karman  Lab 
Main  Stop  301-46 
ATTN:  F.E.C.  Culick 
1201  E.  California  Street 
Pasadena,  CA  91109 

1  California  Institute  of  Technology 

Jet  Propulsion  Laboratory 
ATTN:  LD.  Strand,  MS  512/102 
4800  Oak  Grove  Drive 
Pasadena,  CA  911 09-8099 

1  University  of  Illinois 

Department  of  Mechanical/Industry 
Engineering 
ATTN:  H.  Krier 
144  MEB;  1206  N.  Green  St. 

Urbana,  IL  61801-2978 

1  University  of  Me^sachusetts 

Department  of  Mechanical  Engineering 
ATTN:  K,Jakus 
Amherst.  MA  01002-0014 

1  University  of  Minnesota 

Department  of  Mechanical  Engineering 
ATTN:  E.  Fletcher 
Minneapolis,  MN  55414-3368 


No.  of 

Copies  Organization 

3  Georgia  Institute  of  Technology  School 

of  Aerospace  Engineering 
ATTN:  B.T.  Zim 
E.  Price 
W.C.  Strahle 
Atlanta,  GA  30332 

1  Institute  of  Gas  Technology 

ATTN:  D.  Gidaspow 
3424  S.  State  Street 
Chicago,  IL  60616-3896 

1  Johns  Hopkins  University 

Applied  Physics  Laboratory 
Chemical  Propulsion 
Information  Agency 
ATTN:  T.  Christian 
Johns  Hopkins  Road 
Uurel,  MD  20707-0690 

1  Massachusetts  Institute  of  Technology 

Department  of  Mechanical  Engineering 
ATTN:  T.  Toong 
77  Massachusetts  Avenue 
Cambridge,  MA  02139-4307 

1  Pennsylvania  State  University 

Department  of  Mechanical  Engineering 
ATTN:  V.  Yang 

University  Park,  PA  16802-7501 

1  Pennsylvania  State  University 

Department  of  Mechanical  Engineering 
ATTN:  K.  Kuo 

University  Park,  PA  16802-7501 

1  Pennsylvania  State  University 

Assistant  Professor 
Department  of  Mechanical  Engineering 
ATTN:  Dr.  Stefan  T.  Thynell 
219  Hallowell  Building 
University  Park,  PA  16802-7501 

1  Pennsylvania  State  University 

Director,  Gas  Dynamics  Laboratory 

Department  of  Mechanical  Engineering 
ATTN:  Dr.  Gary  S.  Settles 
303  Mechanical  Engineering  Buldling 
University  Park,  PA  16802-7501 
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1  SRI  International 

Propulsion  Sciences  Division 
ATTN:  Tech  Library 
333  Ravenwood  Avenue 
Menlo  Park,  CA  94025-3493 

1  Rensselaer  Ploytechnic  Institute 
Department  of  Mathematics 
Troy.  NY  12181 

2  Director 

Los  Alamos  Scientific  Lab 
ATTN:  T3/D.  Butler 

M.  Division/B.  Craig 
P.O.  Box  1663 
Los  Alamos,  NM  87544 

1  General  Applied  Sciences  Lab 

ATTN:  J.  Erdos 
77  Raynor  Ave. 

Ronkonkama,  NY  11779-6649 

1  Battelle  PNL 

ATTN:  Mr.  Mark  Gamich 
P.O.  Box  999 
Richland.  WA  99352 

1  Stevens  Institute  of  Technology 

Davidson  Laboratory 
ATTN:  R.  McAlevylll 
Castle  Point  Station 
Hoboken.  NJ  07030-5907 

1  Rutgers  University 

Department  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  S.Temkin 
University  Heights  Campus 
New  Brunswick,  NJ  08903 

1  University  of  Southern  California 

Mechanical  Engineering  Department 
ATTN:  0HE200,  M.  Gerstein 
Los  Angeles.  CA  90089-5199 


1  University  of  Utah 

Department  of  Chemical  Engineering 
ATTN:  A.  Baer 

Salt  Lake  City,  UT  84112-1194 

1  Washington  State  University 

Department  of  Mechanical  Engineering 
ATTN:  C.T.  Crowe 
Pullman.  WA  99163-5201 

1  Alliant  Techsystems,  Inc. 

ATTN:  R.E.  Tompkins 

MN38-3300 

5700  Smetana  Dr. 

Minnetonka.  MN  55343 

1  Alliant  Techsystems,  Inc. 

ATTN:  J.  Kennedy 
7225  Northland  Drive 
Brooklyn  Park,  MN  55428 

1  Science  Applications,  Inc. 

ATTN:  R.B.  Edelman 
23146  Cumorah  Crest  Drive 
Woodland  Hills,  CA  91364-3710 

1  Battelle  Columbus  Laboratories 

ATTN:  Mr.  Victor  Levin 
505  King  Ave. 

Columbus,  OH  43201-2693 

1  Allegheny  Ballistics  Laboratory 

Propulsion  Technology  Department 
Hercules  Aerospace  Company 
ATTN:  Mr.  Thomas  F.  Farabaugh 
P.O.  Box  210 

Rocket  Center.  WV  26726 

1  MBR  Research  Inc. 

ATTN:  Dr.  Moshe  Ben-Reuven 
601  Ewing  St..  Suite  C-22 
Princeton.  NJ  08540 

Aberdeen  Proving  Ground 

1  Cdr,  CSTA 

ATTN:  STECS-PO/R.  HendrIcksen 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  laboratory  urxjertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it 
publishes.  Your  comments/answers  below  will  aid  us  in  our  efforts. 

1 .  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of 
interest  for  which  the  report  will  be  used.)  _ 


2.  How,  specifically,  is  the  report  being  used?  (information  source,  design  data,  procedure, 
source  of  ideas,  etc.)  _ _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or 
dollars  saved,  operating  costs  avoided,  or  efteiencies  achieved,  etc?  If  so,  please 
elaborate.  _ 


4.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports? 
(Indicate  changes  to  organization,  technical  content,  format,  etc.)  _ 


BRL  Report  Number  brl-tr-3362 _  Division  Symbol 

Check  here  if  desire  to  be  removed  from  distribution  list.  _ 

Check  here  for  address  change.  _ 

Current  address:  Organization  _ 

Address  _ 


Department  of  the  army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground.  MO  21005-5066 
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Postage  will  be  paM  by  addressee. 
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U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T 
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